PROCEEDINGS 

OF SCIENCE 



HERWIRI2: CEEX Electroweak Corrections in a 
Hadronic MC 



J2' 



Scott Yost* 1 

The Citadel 

Charleston, SC 29409, USA 
E-mail: |scott . yost@ citadel . edu 



Valerie Halyo * 

Princeton University 
Princeton, NJ 08544, USA 



E-mail: valerieh@princeton . edu 



Miroslav Hejna 

Princeton University 
Princeton, NJ 08544, USA 
E-mail: }nhe jna@princeton . edu 



B.F.L. Ward § 

Baylor University Waco, TX 76798, USA 



E-mail: bfl_ward@baylor.edu 



Reaching the 1 % precision level for W and Z production calculations for the LHC will require a 
mixture of higher order QCD and electroweak corrections. As a first step toward implementing 
the combined QED <_ QCD exponentiation proposed in previous work, we have implemented 
the &{a) electroweak corrections and YFS exponentiation structure of the KK Monte Carlo in 
HERWIG. We discuss the current status of this program and sketch the further developments 
needed to reach the desired precision level. 



10th International Symposium on Radiative Corrections (Applications of Quantum Field Theory to 
Phenomenology ) 
September 26-30, 2011 

Mamallapuram, India BU-HEPP - 11-06 



* Speaker. 

'This work and its presentation were supported in part by D.o.E. grant DE-PS02-09ER09-01 and grants from The 
Citadel Foundation. 

*Work supported in part by D.o.E grant DE-FG02-91ER40671. 
§ Work supported in part by D.o.E. grant DE-FG02-09ER41600. 



<_) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it/ 



HERWIRI2: CEEX Electroweak Corrections in a Hadronic MC 



Scott Yost 



Vector boson production is one of the most important Standard Model processes observed at 
the LHC, and electroweak radiative corrections will be needed for analysis at the percent level. 
Previous studies ^ by some of the authors have found that electroweak corrections alone can 
exceed 1% for some cuts of interest. These studies were based on HORACE, |§] which provides 
state-of-the-art &{a) radiative corrections with a final-state photon shower, and PHOTOS, [^] 
which adds final state photonic radiation. Other programs developed for implementing electroweak 
corrections for hadronic collisions include WINHAC [Q] and ZINHAC[§]. Initial state photonic 
radiation has been considered only in certain MRST partem distribution functions. [^] However, 
none of the most recent PDFs include QED effects. 

For electron-positron colliders, precision electroweak corrections have been implemented in 
the program KKMC 0] using YFS [§] exponentiated multiple-photon radiation for both the initial 
and final state, together with &{a) electroweak corrections in the DIZET []|] package developed 



for ZFITTER[|10p. The DIZET corrections can be applied to any parton-level process, so there is 
no obstacle to extending them to hadronic collisions. YFS exponentiation provides the basis for an 
efficient representation of multi-photon phase space, which can also be implemented for hadronic 
initial and final states. Some of the authors have proposed using a non-abelian extension of YFS 
exponentiation as a basis for an integrated event generator implementing both multi-photon and 
multi-gluon corrections in a unified framework, called QCD&QED exponentiation. JTT| ] 

The collection of programs implementing QCD(g>QED exponentiation has been called HER- 
WIRI, for High Energy Radiation with Infra-Red Improvements, with a version number distin- 
guishing the class of corrections included. The name acknowledges that the initial versions build 



upon the HERWIG [12] parton shower generator. The first to be released, HERWIRI1, [13] imple 



mented IR-improved splitting kernels [14] obtained using the QCD analog of YFS exponentiation. 
This program is publicly available, and tests are in progress. The IR-improved kernels have also 
been implemented [Q in MC@NLO @ and POWHEG[jT7|]. The structure of HERWIRI is not 



tied to a particular shower, and our ultimate goal is a complete shower generator based entirely on 
QCDigiQED exponentiation with exact @(a%,a s a, a 2 ) residuals. [ ]T8| ] 



The second version, HERWIRI2, implements the electroweak radiative corrections of KKMC 
in a hadronic shower generator. This note describes a version of HERWIRI2 which is presently 
nearing completion. 1 HERWIRI2 is motivated by the successful application of YFS exponentiation 
in BHLUMI[|l9|], BHWIDE@, KoralZ[|I|], KoralW[j22|], YFSWW3[||], KKMC, and related 



programs for LEP physics. All of these programs benefit from a very efficient representation of N- 
photon phase space, with complete control over the soft and collinear singularities for an arbitrary 
number of photons. Real and virtual IR singularities cancel exactly to all orders. The non-abelian 
extension to QCD® QED should have similar advantages. 

While based on HERWIG, HERWIRI2 is largely independent of the underlying shower. HER- 
WIG generates the parton momenta and shower, and HERWIRI2 passes the generated hard process 
momenta to KKMC routines to add photons and electroweak corrections. Although KKMC was 
developed for e + e~ collisions, it was designed to be extended to more general processes, so the 
ability to select quarks as the incoming state already exists in all but the lowest-level generation 
routines. 



HERWIRI2 does not incorporate HERWIRI1, although the two programs can be used in combination. 
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KKWIC is a precision generator for e + e~ — > ff + ny, f = n,r,d,u,s,c,b for CMS energies 
from 2m t to 1 TeV. The precision tag for LEP2 was 0.2%. ISR and FSR y emission is calculated 
up to iff (a 2 ), including interference. The MC structure is based on YFS exponentiation, including 
residuals calculated perturbatively to the relevant orders in a k L l . (L = \n(s/m 2 )). Exact collinear 



bremsstrahlung is implemented for up to three photons. Electroweak corrections [ ]24| , [25| , [26[ ] are 
included via DIZET 6.2 and beamstrahlung can be modeled over a wide range of energies via 
a built-in or user-defined distribution. 

There are two modes of operation: exclusive exponentiation (EEX) and coherent exclusive ex- 
ponentiation (CEEX). EEX applies YFS exponentiation to differential cross-sections, while CEEX 
applies it at the amplitude level. The CEEX mode is written in a manner that is most readily 
extended to quark scattering, so it is taken as the basis for HERWIRI2. The orders of residuals 
included in CEEX mode are a, aL, a 2 L 2 , and a 2 L. 

CEEX was introduced for pragmatic reasons, because the traditional (EEX) exponentiation 
of spin-summed differential cross sections suffered from a proliferation of interference terms in 
processes with multiple diagrams, limiting its utility in reaching the desired 0.2% precision tag 
for LEP2.[|J CEEX works at the level of spinor helicity amplitudes, greatly facilitating the cal- 
culation of effects such as ISR-FSR interference, which are included in KKMC, and therefore in 
HERWIRI2. CEEX is maximally inclusive: all real photons radiated are kept in the event record, 
no matter how soft or collinear. There is no need to integrate out a region of soft phase space, 
because the exponentiated amplitudes are well-behaved at k = 0. 

The CEEX cross section for qq — > f f has the form [[7|] 
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The YFS form factor is 
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where Qj, Qf are initial and final parton charges, and 

Y(pi,pj,E min ) = 2aB(pi,pj,E^ m ) + 2aReB(pi,pj), 
with real and virtual form factors defined respectively by 
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The «-photon helicity-spinor amplitude can be expanded in terms of order a r having the form 
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with residual spinor amplitudes J3^ and complex soft photon factors S^j with the property 
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where Q,p a ,Pb belong to the initial or final fermions depending on the partition 

KKMQ, incorporates the DIZET library (version 6.21) from the program ZFITTER. [1C] The 
7 and Z propagators are multiplied by vacuum polarization factors: 
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Vertex corrections are incorporated into the coupling of Z to / via form factors in the vector cou- 
pling: 
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Box diagrams contain these plus a new angle-dependent form-factor in the doubly-vector compo- 
nent: 

(zj) (zj) _ 7f ?f) _ 2T^Q fF y\s) - 2Q t T^%) +4Q iQf F^ (s,t) 
8v 8v ~ sh?(2M ' 

The correction factors are calculated at the beginning of a run and stored in tables. 

The Drell-Yan cross section with multiple-photon emission can be expressed as an integral 
over the parton-level process qi{p\)q~i{p2) — > f{pi)f{pA) + ny(k), integrated over phase space and 
summed over photons. The parton momenta p\ , pi_ are generated using parton distribution functions 
giving a process at CMS energy q and momentum fractions x\ ,X2 such that q 2 = x\x 2 s\ 



0dy 



dx\ dx2 
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(12) 



where the final state phase space includes p^,p4 and i = 1, • • • ,n and multiple gluon radiation + 
hadronization is included through a shower. 

The parton-level cross section Oj(q 2 ) can be calculated by KKMC, which integrates over a 
final space phase space with two fermions and an arbitrary number of photons: 



Oi(q 2 ) = £ / dPS 2+ nOi(p,k) 

n=0 J 



(13) 



HERWIRI2 uses HERWIG 6.5 as the shower generator, which creates the hard process first, 
at Born level, in subroutine HWEPRO (HWHDYP), and then passes it to the cascade generator 
HWBGEN. HERWIRI2 finds the Z/y* and the partons interacting with it in the event record. The 
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initial partons define p\, p2, which are transformed to the CM frame and projected on-shell to 
create a starting point for KKMC, which generates the final fermion momenta p$,p4 and photons 
ki (both ISR and FSR.) The generated particles are transformed back to the lab frame and placed in 
the event record. 

In addition to the basic DY process, HERWIG generates "Compton" events g + q — >■ q + Z/j* ■ 
About 10% of the events have this form. This is factorized into gluon emission times a hard 
EW process at a shifted value of q 2 . These have a different profile in the event record, but can be 
processed by KKMC as well. There is also a third class of events with the emission of an additional 
hard gluon. About 1% of the events have this form, and also have a significant shift of the Z energy 
from its generation scale. 

With a change of variables, the Drell-Yan cross section in HERWIG can be expressed as 

Ody = [— — y,fi(q,xi)f 1 (q,x 2 )o i (q 2 )8(q 2 -x l X2s) 
J X\ X2 *Y 

= r™dqP(q) \\ W®, (q 2 , Xl ) = (W HW ) (14) 

where P(q) is a normalized, integrable, crude probability distribution for q, Pi is the crude proba- 
bility of generating parton i, and Whw is the HERWIG event weight. This weight depends only on 
the hard Born cross section and is not altered by the shower. 
The crude probability distributions used by HERWIG are 

The HERWIG event weight is 

Whw = I< w , < w = j^f^xiMq^ln (^J oj^(q 2 ) (16) 
and the corresponding probability for selecting parton i is 

Pi = </WW (17) 

We have chosen to introduce electroweak corrections in a minimally invasive way, incorporat- 
ing them in a form factor 

KKMC will calculate the EW form factor, and multiply it by the HERWIG Born cross section. 

(jW ( a 2) 

(Thw+ew = (Wrot) , W Tot = F^(q 2 )W B w = W B w ^ ' ■ (19) 
The KKMC cross section is calculated using a primary distribution 
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with 



2a - 
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to generate the factor v giving the fraction of s remaining after ISR photon emission, sx = s(l — v). 
The KKMC cross section is 

a(q 2 ) = f doJ-p^ d -p^ = a Pn (W Cm W Mod ) . (22) 

Wcm is calculated during ISR generation and WMod is generated after sx is available. 

The HERWIG and KKMC weights are combined to calculate the total HERWIRI2 weight, 
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This average com/g? be calculated using a joint probability distribution for q and v, D(q,v) 



P(q)da Pli /dv, withP(<7) from HERWIG. An adaptive MC (S. Jadach's FOAM @|) could calculate 
the normalization of the distribution at the beginning of the run, in a similar manner to how KKMC 
presently integrates the one-dimensional primary distribution. To account for beamsstrahlung, 
KKMC already permits such a user-defined distribution, in up to three variables. However, as a 
first step, we have tried to run HERWIRI2 using KKMC's one-dimensional primary distribution. 
This requires fixing an overall scale qo to initialize KKMC (e.g., qo = M-£). 

The built-in primary distribution for electrons at scale qo can be used for the low-level genera- 
tion of v. The transformation from this distribution to a distribution at HERWIG's generated scale 
q for quark i is then obtained by a change of variables: 

with 

^tf,v)_^0<^.tf(l-v» | (25) 



where 

Wy ■■£-[ — ) (26) 



mm 
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The / factors are calculated using q 2 /mf for parton i and ql/m 2 for the electron. 

Shuffling the numerators and denominators about gives the expression used in HERWIRI2: 



a To t = (WuwW M odW Karl W w W Y ) (27) 



with two new weights 
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HERWIRI2 is still under development, so any numerical results must be treated as prelimi- 
nary. A run for pp collisions at 10 TeV with the Z/y* invariant mass bounded by 40 GeV and 140 
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GeV, using HERWIG 6.520 default parameters and PDFs, yields a cross-section of 1183.7 ± 1.3 
pb, compared to 1098.8 ± 1.0 pb for HERWIG alone, giving an electroweak contribution of 7.7%. 
Turning on ISR gives a much wider weight distribution and consequently, greatly reduced effi- 
ciency. Preliminary results give a cross-section of 1212 ± 109 pb, showing an additional 2.4% 
contribution from ISR. 

Work is in progress to optimize MC generation in the presence of ISR. Once HERWIRI2 
is complete and thoroughly tested, it will be compared to other available hadronic/electroweak 
generators. In particular, it will be interesting to see the effect of initial state radiation, which is not 
present in the other programs, but appears to enter at the 2-3% level, making it crucial to precision 
calculations. 
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